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The mixed spin-1/2 and spin-^* {S > 1/2) Ising model on a rope ladder 
is examined by combining two exact analytical methods. By the decoration- 
iteration mapping transformation, this mixed-spin system is firstly trans- 
formed to a simple spin-1/2 Ising model on the two-leg ladder, which is 
then exactly solved by the standard transfer-matrix method. The thermal 
variations of the zero-field specific heat are discussed in particular. 



PACS numbers: 05.50.-Hq, 75.40.Cx 

1. Introduction 

One-dimensional Ising models are of fundamental scientific interest partly 
on behalf of their exact solubility [T] and partly due to the fact that they 
represent the simplest lattice-statistical models, which have found rich appli- 
cations in seemingly diverse research areas [2]. In the present work, we shall 
provide the exact solution for the mixed spin-1/2 and spin-S {S > 1/2) Ising 
model on a rope ladder, which represents a rather novel magnetic structure 
emerging in several insulating bimetallic coordination compounds [3j. 

2. Model and its exact solution 

Consider the spin-1/2 Ising ladder of 2N sites with each of its horizontal 
and vertical bonds occupied by the decorating spin S {S > 1/2), as it is 
schematically shown in Fig. [TJ For further convenience, the Hamiltonian of 
the model under investigation can be written as a sum of site Hamiltonians 

3N 

Ti = ^Hk, 'Hk = -JSk{o-ki + (Tk2) - DSl, (1) 

k=l 
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Fig. 1. A small fragment of the mixed-spin Ising model on the rope ladder. The 
full circles denote lattice positions of the vertex spins Uj = ±1, while the empty 
ones label lattice positions of the decorating spins Sk ~ ^S, —S + 1, . . . , +S. 



where each site Hamiltonian 7ik involves all the interaction terms of one 
decorating spin. Above, cr^^ = ±1 (7 = 1,2) denote two vertex spins of 
the rope ladder that are nearest neighbours to the decorating spin Sk = 
—S, —S + 1, . . . , +S. The parameter J labels pairwise exchange interaction 
between the nearest-neighbour spin pairs, while the parameter D measures 
a strength of the single- ion anisotropy acting on the decorating spins. 

The crucial step of our procedure lies in the calculation of the partition 
function, which can be firstly partially factorized to the expression 

3N +S 

^ = E n E c°sh [pjn (cTfci + ak2)] , (2) 

{a-j} j=l n=-S 

where /3 = l/(kBT), k-Q is Boltzmann's constant, T is the absolute temper- 
ature and the symbol J2{cr} denotes the summation over all possible spin 
configurations of the vertex spins aj. Next, it is advisable to employ the 
generalized decoration-iteration mapping transformation 

+s 

^ exp (^/3Dn'^'j cosh [/3Jn {aki + cjfc2)] = ^exp {l3RakiCFk2) , (3) 

n=-S 

which physically corresponds to removing all the interaction parameters 
associated with the decorating spins and replacing them by a new effective 
interaction R between the remaining vertex spins. Since there are four 
possible spin states for each pair of vertex spins aki and ak2 and only two of 
them provide independent equations from the transformation formula ([3]), 
the mapping parameters A and R have to meet the following conditions 

A^=l ^ exp(/3Dn2)cosh(2/?Jn) I ( ^ exp(/3L»n2) | , (4) 

\n=-5 / \n=-5 / 

(5R = -\ni exp{(3Dn^) cosh(2/3Jn) | - - In [ ^ exp{PDn'^)\. (5) 

^ \n=-S / ^ \n=-S / 
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By substituting the transformation ([3]) into Eq. ([2]) one gains the relation 

2(/3,J,i?) =A3^2o(/3,i?), (6) 

which connects the partition function Z of the mixed-spin Ising model on 
the rope ladder with the partition function Zq of the simple spin-1/2 Ising 
model on two-leg ladder with the effective nearest-neighbour interaction R. 
It is worthwhile to remark that the partition function of the latter model 
can be rather easily calculated by applying the transfer-matrix method [1], 
which yields in the thermodynamic limit the following exact result 

Zo = |cosh(3^ii) + cosh(/3ii) + \l [sinh(3/?i?) - s\nh[[3R)f + 4| . (7) 

Bearing this in mind, the exact solution for the partition function of the 
mixed-spin Ising model on the rope ladder is formally completed as it can 
be simply achieved by a mere substitution of the corresponding partition 
function ([7]) to the relationship ([6]), whereas both mapping parameters A 
and R have to be taken from Eqs. (jlj) and ([5]), respectively. 

3. Results and discussion 

Now, let us focus on thermal variations of the zero-field specific heat with 
the aim to shed light on how this dependence changes with the single-ion 
anisotropy. Before discussing the most interesting results, it should be men- 
tioned that the zero-field specific heat can readily be obtained from the exact 
result for the partition function ^ with the help of basic thermodynamical- 
statistical relations. Even although the final expression is too cumbersome 
to write it down here explicitly, it is worthy to note nevertheless that it 
implies an independence on a sign of the exchange interaction J. In this re- 
spect, the temperature dependences of the zero-field specific heat discussed 
below remain in force regardless of whether ferromagnetic (J > 0) or ferri- 
magnetic ( J < 0) spin system is considered. 

For illustration. Fig. [2^) and b) depict thermal variations of the zero- 
field specific heat for the rope ladder with the decorating spins S = 1 and 
S = 3/2. As one can see, the same general trends can be observed in both 
figures: the round maximum diminishes and shifts towards lower tempera- 
tures by decreasing the single- ion anisotropy. However, the most interesting 
temperature dependences can be found at sufficiently negative single-ion 
anisotropies close to D/ J = —2.0 for the case with S = 1 and respectively, 
D/J = —1.0 for the case with S = 3/2. In the vicinity of both these 
boundary values, the double-peak specific heat curves (low-temperature 
peaks are for clarity shown in the insets) appear due to the competition 
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Fig. 2. The temperature dependences of zero-field specific heat for several single- ion 
anisotropics and two different spin values S* = f (Fig. [5^) and S ~ 3/2 (Fig.Os). 



between two different spin configurations sufficiently close in energy. Actu- 
ally, it can easily be verified that the general condition ensuring an ener- 
getic equivalence between two different spin states of decorating spins reads 
Ds^S-i = —2J/{2S — 1). It should be also mentioned that the zero-field 
specific heat curves with a single maximum are recovered upon further de- 
crease of the single- ion anisotropy. In conclusion, it seems valuable to extend 
our calculation to more general case with the non-zero magnetic field. 
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